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We review the equation of state for dark energy in modified gravity theories. In 
particular, we summarize the generic feature of the phantom divide crossing in the past 
and future in viable f(R) gravity models. 
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To understand the late time acceleration universe one of the interesting 
possibilities is to consider a modified gravitational theory, such as f(R) grav- 
ity To build up a viable f(R) gravity model, one needs to satisfy the fol- 
lowing conditions: (a) positivity of the effective gravitational coupling, (b) sta- 
bility of cosmological perturbations "I, (c) asymptotic behavior to the standard 
A-Cold-Dark-Matter (ACDM) model in the large curvature regime, (d) stabil- 
ity of the late-time de Sitter point ^1, (e) constraints from the equivalence prin- 
ciple, and (f) solar-system constraints Several viable mode ls have been con- 
structed in the literature, such as the following popular onesP^ ^^^^^^ I 

model f(R) Constant parameters 

(i) Hu-Sawicki R - C ^W*$ S X ci, c 2 , p(> 0),i? HS (> 0) 



(ii) Starobinsky R + XR§ 



1 



A(> 0), n(> 0), Rs 



(hi) Tsujikawa R - (iRr tanh (-J^ fi(> 0), R T {> 0) 

(iv) Exponential R - f3R E (l - e - R / RE ) f3, R E 

Recently, the cosmological observational data^ seems to indicate the crossing of 
the phantom divide ?«de = — 1 of the equation of state for dark energy in the near 
past. To understand such a crossing, many attempts have been made. The most 
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noticeable one is to use a phantom field with a negative kinetic energy term H^l 
Clearly, it surfers a serious problem as it is not stable at the quantum level. On the 
other hand, the cr ossing of the p hantom divide can also be realized in the above 
viable f(R) models ' ' l | l^ | io | lb| without violating any stability conditions. This is 
probably the most peculiar character of the modified gravitational models. Oth er 
f{R) gravity models with realizing a crossing as well as multiple crossings of 
the phantom boundary have also been examined. 

In this talk, we would like to review equation of state in /(i?) gravity. In par- 
ticular, we show that the viable f(R) model s gener ally exhibit the crossings of the 
phantom divide in the past as well as future I 12 | 13 | 

The action of f(R) gravity with matter is given by 

I = J ^ X ^^^~2~T + Imatter(9nu, Y matter) ) 00 

where g is the determinant of the metric tensor g^ u , /matter is the action of matter 
which is assumed to be minimally coupled to gravity, i.e., the action / is written 
in the Jordan frame, and T ma tter denotes matter fields. Here, we use the standard 
metric formalism. By taking the variation of the action in Eq. (J]} with respect to 
g^ v , one obtains^ 

FG^ = ^matter) _ ^ {FR _ f) + _ g ^ up ; (2) 

where G M „ = R^ w — (1/2) g^R is the Einstein tensor, F(R) = df(R)/dR, V M is the 
covariant derivative operator associated with g^ v , □ = g^ u V ^V„ is the covariant 
d'Alembertian for a scalar field, and y^™ attor ) j s the contribution to the energy- 
momentum tensor from all perfect fluids of matter. 

From Eq. ([2]), we obtain the following gravitational field equations: 

3FH 2 = K, 2 p M + X - (FR - f) - 3HF , 

-2FH = k 2 (pm +P M ) + F-HF, (3) 

where H = d/a is the Hubble parameter, the dot denotes the time derivative of 
d/dt, and pm and Pm are the energy density and pressure of all perfect fluids of 
matter, respectively. 

The equation of state for dark energy is given by 



w DE = Pde/pde, (4) 



where 



Pde — 9 



Pde — r 

K 



- {FR - f) - 3HF + 3(1-F)H 2 
-- (FR - f)+F + 2HF - (1 - F) (2H + 3H 2 



(5) 
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In Figs. 1, 2 and 3, we depict the evolution of wde, future evolutions of 1 + wde, 
and H = H — Hf with H = H/Hq and Hf = H(z = —1)/Hq, as functions of the 
redshift z = 1/a — 1 in (i) Hu-Sawicki model for p = 1, c\ = 2 and C2 = 1, (h) 
Starobinsky model for n = 2 and A = 1.5, (hi) Tsujikawa model for /i = 1 and 
(iv) the exponential gravity model for j3 = 1.8, respectively, and the subscript 'f 
denotes the value at the hnal stage z = — 1. Note that the present time is z — 
and the future is — 1 < z < 0. The parameters used for each model in Figs. 1-3 
are the viable ones EE12QI. Several remarks are as follows: (a) the qualitative results 
do not strongly depend on the valu es o f the parameters in each model; (b) we have 
studied the Appleby-Battye model ^H, which is also a viable f{R) model, and we 
have found that the numerical results are similar to those in the Starobinsky model 
of (ii) as expected. 

We note that the present values of wde(z = 0) are -0.92, -0.97, -0.92 and -0.93 




-1 0123456 -1 012345678 




Fig. 1. Evolutions of the equation of state idde &s functions of the redshift z in (i) Hu-Sawicki 
model for p = 1, ci = 2 and C2 = 1, (ii) Starobinsky model for n = 2 and A = 1.5, (iii) Tsujikawa 
model for fi = 1 and (iv) the exponential gravity model for /3 = 1.8, respectively, where the thin 
solid lines show iode = ~ 1 (cosmological constant). 
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for the models of (i)-(iv), respectively. These values satisfy the present observa- 
tional constraints k*2l Moreover, a dimensionless quantity H 2 / /3) can be 

determined through the numerical calculations, where pS' is the energy density 
of non-relativistic matter at the present time. If we use the observational data on 
the current density parameter of non-relativistic matter Q$ = Pm / P^ lt — 0.26 
with p^) it = 3Hq I k 2 ^21, we find that the present value of the Hubble parameter 
H = H(z = 0) is H = 71km/s/Mpc 1221 for all the models of (i)-(iv). Further- 
more, Hf = 0.80 , 0.85 , 0.78 and 0.81 for the models of (i)-(iv), respectively, where 
Hi = H(z = -l)/H . 

It is clear from Figs. 1-3 that in the future (—1 < z < —0.74), the cross- 
ings of the phantom divide are the generic feature for all the existing viable f(R) 
models. By writing the first future crossing of the phantom divide and the first 
sign change of H from negative to positive as z = z cmsa and z = z p , respectively, 
we find that (z cross ,z p ) a = (-0.76,-0.82),, (-0.83, -0.98)«, (-0.79, -0.80)^ and 
(—0.74, — 0.80)j„, where the subscript a represents the ath viable model. The values 
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of the ratio 3 = f2 m /f2DE at z = z cmss and z = z p are (S(z = z C ross), 3(z = z p )) a — 
(5.2 x 10" 3 , 2.1 x 10~ 3 );, (1.7 x 10~ 3 ,4.8 x 10"%, (4.1 x 10~ 3 ,3.1 x lO" 3 )™ and 
(6.2 x 10~ 3 ,2.8 x 10~ 3 )„, where ft DE = Pue/p^X and = Pm/ Pelt are tne den- 
sity parameters of dark energy and non-relativistic matter (cold dark matter and 
baryon), respectively. As z decreases (—1 < z < —0.90), dark energy becomes much 
more dominant over non-relativistic matter (S = f2 m /f2DE ^ 10~ 5 ). As a result, one 
has wde ~ w e s = —1 — 2H / (3i/ 2 ) = -Ptot/ptot, where w e s is the effective equation 
of state for the universe, and ptot = Pde + Pm + Pi and Ptot = -Pde + Pr are the total 
energy density and pressure of the universe, respectively. Here, p m ( r ) and P r are the 
energy density of non-relativistic matter (radiation) and the pressure of radiation, 
respectively. The physical reason why the crossing of the phantom divide appears 
in the farther future (— 1 < z < —0.90) is that the sign of H changes from negative 
to positive due to the dominance of dark energy over non-relativistic matter. As 
iude ~ i" c n in the farther future, u>de oscillates around the phantom divide line 
wde = — 1 because the sign of H changes and consequently multiple crossings can 
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be realized. 

Finally, we mention that in our numerical calculations, we have taken the initial 
conditions of zq = 8.0, 8.0, 3.0 and 3.5 for the models of (i)-(iv) at z = zq, respec- 
tively, so that RF'(z — zq) ~ 10~ 13 with F 1 = dF/dR, to ensure that they can be 
all close enough to the ACDM model with RF 1 = 0. 

In this talk, we have explored the past and future evolutions of u>de in the viable 
f{R) gravity models and explicitly shown that the crossings of the phantom divide 
are the generic feature in these models. We have demonstrated that in the future 
the sign of H changes from negative to positive due to the dominance of dark energy 
over non-relativistic matter. This is a common physical phenomena to the existing 
viable f(R) models and thus it is one of the peculiar properties of f(R) gravity 
models characterizing the deviation from the ACDM model. 
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